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ABSTRACT 

We report the discovery of HAT-P-56b by the HATNet survey, an inflated hot Jupiter transiting a 
bright F type star in Field 0 of NASA’s K2 mission. We combine ground-based discovery and follow-up 
light curves with high precision photometry from K2, as well as ground-based radial velocities from 
TRES on the FLWO 1.5 m telescope to determine the physical properties of this system. HAT-P-56b 
has a mass of 2.18 Mj, radius of 1.47 i?j, and transits its host star on a near-grazing orbit with a period 
of 2.7908 d. The radius of HAT-P-56b is among the largest known for a planet with Mp > 2Mj. The 
host star has a U-band magnitude of 10.9, mass of 1.30 Mq, and radius of 1.43i?©. The periodogram 
of the K2 light curve suggests the star is a 7 Dor variable. HAT-P-56b is an example of a ground-based 
discovery of a transiting planet, where space-based observations greatly improve the confidence in the 
confirmation of its planetary nature, and also improve the accuracy of the planetary parameters. 

Subject headings: planetary systems — stars: individual (HAT-P-56) — techniques: spectroscopic, 
photometric 


1. INTRODUCTION 

Searching for transits is one of the most productive 
methods of detecting planets outside of our solar sys¬ 
tem. Transiting exoplanets (TEPs) around bright stars 
are of particular interest as they are the best targets 
for follow-up observations to investigate their detailed 
orbital geometries, atmospheric properties and chemical 
compositions. 

Recent years have seen a rapid increase in the dis¬ 
covery rate of transit exoplanet s (TEPs). The Ke¬ 
pler Mission (jBorucki et al.l[20?)^ has identified ^4000 
planetary candidate s in its four years of operations 
(jMullallv et al.l[2015f) . Approximately a hundred of the 
Kepler planets hav e their mass measu r ed using stel¬ 
lar radial velocities (iBatalha et al.l (120111) , iGautier et al.l 
(I2012D , IMarcv et al.l ( 20141) . and many oth ers), or numer¬ 
ically modeled trans i t tim i ng variations (iLissauer et HI 
(|20HD . [Carter et HI (l2012D . fWu &: LithwickI ( 20131 ) and 
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many others). The majority of Kepler planet candidate 
hosting stars are, however, often faint, complicating both 
the accurate planetary mass measurements and also the 
study of their physical nature. For example, to date, 
the majority of the planets with precise mass measure¬ 
ments are still from the ground-based wi de field transit¬ 
ing p lanet surveys such as S uper WASP (iPollacco et ^ 
HATNet/HATSout h (IBakos et al.l 12004 I20I3D 
KELT (jPepper et al.ll200^ . and others. 

Traditionally, confirmation of planets from ground- 
based transit surveys relies on high precision photometric 
follow-up from one or two meter class ground-based tele¬ 
scopes. These follow-up observations ensure the robust¬ 
ness of the detections, help to rule out false positives due 
to binaries, and further constrain the planet parameters. 
However, due to the constraints of ground-based obser¬ 
vations, it’s usually difficult to observe the full orbital 
phase of a transiting planet candidate, including primary 
and secondary transits. The presence or lack of deep sec¬ 
ondary eclipses, or out-of-transit variations is often im¬ 
portant information for ruling out blended stellar eclips¬ 
ing binary scenarios. The quality of ground-based pho¬ 
tometry may also be reduced due to poor weather con¬ 
ditions, site and target restrictions, and airmass trends. 

Due to the loss of two reaction wheels during the main 
mission, the Kepler spacecraft entered a new observation 
phase. This successor miss ion of Kepler, called the K2 
mission (iHowell et al.ll201^ . covers a much larger area of 
sky with a step and stare strategy. The Guest Observa¬ 
tion mode of this mission provides a great tool to follow 
up the ground-based planetary candidates, and enables a 
new method of discovering TEPs through a synergy be¬ 
tween ground- and space-based transit surveys. K2 ob¬ 
servations allow us to obtain high precision light curves 
of candidates over a continuous and relatively long time- 
baseline covering many transit events due to a planet. 

Once the planet is confirmed, having high precision 








































2 


Huang et al. 


light curves from space observations will also enable fur¬ 
ther characterization of the planet. Conhrmed planets 
from the ground have bee n recognized as va l uable targets 
to be followed up by K2 (jBakos et al.ll2ni5l : iBrown et al.l 
l2014fl . The planets discovered by ground based TEP 
surveys, and also obse rved by the original K epler Mis- 
sio n, such as TrEs -2b (|0’Donovan et ai1l2006f). HAT-P - 
7b (IPal et al.ll2008ll . and HAT-P-llb (iBakos et^lWfill . 
are among the best studied planets. With the high 
precision Kepler light curves, many effects due to the 
planets can be measured, such as the orbital phase 
varia ti on and occultation o f the planet (iBarc lav et all 
2012t Dackson et al.l 120121: iKipping fc BakosI l2011bl : 
Welsh et al.ll2010t iBorucki et al.l 1200911 . The spin-orbit 
obliquity angle of HAT-P-llb has also been constrained 
using the sta rspot crossing events observed in the Kepler 
light curves (jSanchis-Oieda fc Winnl[20Tll : iDeming et al.l 
l2011fl . Detailed modeling of the transit shape using the 
Short Cadence photometry of Kepler can also be used 
to determine the probable gravity da rkening effect, and 
constrain the oblateness of the planet (iMorris et al.ll20l5 
lMasudall2015l: f^u et al.ll2014l l. 

Following up candidates from ground-based transit 
surveys is an efficient way to utilize the K2 observational 
resources. Due to constraints from the new mode of op¬ 
eration, the number of targets observed by K2 per held is 
much smaller than for the original Kepler held. The high 
priority planetary candidates from ground-based surveys 
are pre-selected targets which are known to show tran¬ 
sits, and which have already been vetted against various 
false positive scenarios. Selecting these candidates to fall 
on K2 “postage stamps” is one way to increase the con¬ 
hrmed planet yield from this mission. 

In this paper we present the discovery (see Figure[T|) of 
a transiting planet, HAT-P-56b, in the K2 Campaign 0 
held. This planet was originally identihed as a HATNet 
(jBakos et al.ll2004ll planetary candidate, was followed up 
by the TRES spectrograph on the FLWO 1.5 m telescope, 
and also by the KeplerCam imager on the FLWO 1.2 m 
telescope. Encouraged by these initial results, all point¬ 
ing toward a bona fide planet orbiting the host star HAT- 
P-56, the target was proposed for K2 observations. In¬ 
deed, the very high quality photometric observations of 
K2 conhrmed the transit, and also eliminated most of the 
possible blend scenarios. We then continued following-up 
HAT-P-56 with the TRES spectrograph, so as to deter¬ 
mine the mass of the orbiting body. 

In ^ we summarize the detection of the photometric 
transit signal in the HATNet light curve, follow-up pho¬ 
tometry from the ground and from K2 campaign 0, and 
our spectroscopic follow-up. Analyses of the results are 
presented in fjS] We show in jjUthat HAT-P-56b is one 
of the most inhated objects observed that belong to the 
massive hot Jupiter population (M > 2Mj), even when 
we take into account of the amount of irradiation from 
the host star. 


2. OBSERVATIONS 

2 .1. Ground Based Photometry 

All time-series photometric data that we collected for 
HAT-P-56 are provided in Table [TJ We discuss these ob¬ 
servations below. All of our discovery data and follow-up 
data are publicly available in electronic format as ma- 
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Fig. 1. — The sequence of observations that lead to the discovery 
and confirmation flow of HAT-P-56b. 

chine readable tables, and also via HATNet website0. 

2.1.1. Photometric detection 

The star HAT-P-56 was observ ed by the HATNet w ide- 
held photometric instruments (IBakos et al.l 120041) be¬ 
tween the nights of UT 2011 October 14 and UT 2012 
May 3. A total of 6509 observations of a 10?6 x 10?6 
held centered at R.A. = 06*'‘'24™''', Dec. = -1-30° were 
made with the HAT-6 telescope in Arizona, and 4194 
observations of this same held were made with the HAT- 
9 telescope in Hawaii (the count is after hltering some 40 
outlier measurements). We used a Sloan r hlter an d an 
exposure time of 180 s. Fohowing IBakos et al.l (|2010ll and 
iKovacs et al.l (1200511 . we reduced the images to trend- 
hltered light curves for the ^ 124,000 stars in the Held 
with r < 14.5 mag, achieving a point-to-point r.m.s. pre¬ 
cision of 3.5mmag for the brightest non-saturated stars 
with r ^ 10 mag (for HAT-P-56 the r.m.s. of the resid¬ 
uals from our best-ht transit model is 6.7mmag). We 
searched these light curves for periodic tra nsit signals us¬ 
ing th e Box Least Squares algorithm IBLS: IKovacs et all 

I200I . 

A total of 29 candidate transiting planets were identi¬ 
hed from these light curves, including HAT-P-56, whose 
phase-folded HATNet light curve we show in Figure [2J 
Reconnaissance spectroscopy and photometric follow-up 
observations have been carried out for most of these can¬ 
didates, based on which we have rejected 21 of them 
as false positives. In addition to HAT-P-56, which we 
conhrm as a planetary system in this paper, one other 
planet has been c onhrmed from this held (HAT-P-54b; 
IBakos et al.ll20igi . Six of the candidates remain active. 

2.1.2. Photometric follow-up with KeplerCam 

Initial photometric follow-up observations of HAT- 
P-56 were carried out with KeplerCam on the Fred 
Lawrence Whipple Observatory (FLWO) 1.2 m telescope. 
We observed a single transit ingress on the night of UT 
2013 March 25 using an Aband hlter and an exposure 
time of 10 s. The images w ere reduced to a light curve fol¬ 
lowing [Bakoi^F^I] (1201011 . including external parameter 
decorrelation performed simultaneously with the transit 
ht to remove systematic trends; this trend-corrected light 
curve is shown in Figure [S] The r.m.s. of the residuals 
from our best-ht model is 2.2mmag for these data. 

® http://hatnet.org 
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Fig. 2.— HATNet light curve of HAT-P-56 phase folded with 
the transit period. The top panel shows the unbinned light curve, 
while the bottom shows the region zoomed-in on the transit, with 
dark filled circles for the light curve binned in phase with a binsize 
of 0.002. The solid line shows the model fit to the light curve. 

2 .2. Spectroscopy 

We carried out spectroscopic observations of HAT-P- 
56 between UT 2012 October 31 and UT 2014 November 
25 using the Ti l lingh ast Reflector Echelle Spectrograph 
(TRES; lEi'iresTl [200l on the 1.5 m Tillinghast Reflec¬ 
tor at ELWO. We extracted spectra from the images and 
measured initial RVs, bisector s pans (BSs) and stellar a t- 
mospheric parameters following iBuchhave et al.l (l2010f) . 

The first three TRES observations made between UT 
2012 October 31 and UT 2012 November 7 using a short 
exposure time of 360 s resulted in a S/N per resolution 
element of 38-45, and were used for reconnaissance pur¬ 
poses. Based on these three reconnaissance spectra we 
found that HAT-P-56 is an E dwarf star with a fairly 
rapid projected rotation velocity of ^ 40kms“^, and no 
evidence for additional stellar components in the spec¬ 
trum. The three RV measurements were consistent with 
no variation above 600ms“^ (3 (t upper limit). A single 
APO 3.5m/ARCES spectrum of HAT-P-56 was also ob¬ 
tained for reconnaissance on UT 2012 November 7. A 
160 s exposure time was used to achieve a S/N per reso¬ 
lution element of 28.6. We reduced the observation to a 
wavelength-calibrated spectrum with the IRAP echelle 
packagcj and used the Stellar Param eter Classification 
program ('SPCi lBiichhave et al.ll201^ to measure the at¬ 
mospheric parameters and radial velocity from the spec¬ 
trum. This spectrum was also found to be single-lined, 
had an RV consistent with the TRES measurements, 
and also indicated that the target was a F dwarf with 
usinf = 40kms“^. 

Following the K2 observations (Section I2.3|l . which 
showed clean transits with no evidence of secondary 
eclipses or strong out-of-transit variability in phase with 
the transits, we resumed spectroscopic monitoring with 
TRES on the FLWO 1.5 m telescope, now with the aim of 
measuring the mass of the planet by detecting the orbital 
motion of the host star. A total of 18 high S/N (ranging 
from 50-113) observations were collected between UT 
2014 October 2 and UT 2014 November 15. We mea¬ 
sured stellar atmospheric parameters from these spectra 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3.— Unbinned transit light curve for HAT-P-56, acquired 
with KeplerCam at the FLWO 1.2 m telescope on the night of 2013 
March 25. The light curve was corrected for trends, which were 
fit simultaneously with the transit model. Our best fit from the 
global modeling described in Section 0 is shown by the solid line. 
Residuals from the fit are displayed below the original light curve. 
The error bars represent the photon and background shot noise, 
plus the readout noise. 

using SPC, a nd carried out a mul ti-order velocity anal¬ 
ysis following iBiervla et al.l ()2014f) to measure the RVs 
relative to one of the observed spectra and spectral line 
bisector spans (BSs). Table [2] provides the measurements 
extracted from these spectra, where we exclude the re¬ 
connaissance spectra that were deemed to be of too low 
S/N to be used for our final characterization of the orbit 
and the stellar atmospheric properties. The phase-folded 
RVs and BSs are shown in Figure |4] together with our 
best-ht circular and eccentric orbit models. We show in 
Section [321 that the above RV measurements confirm the 
planet nature of HAT-P-56b. 

2.3. K2 photometry follow up 

Encouraged by the HATNet, FLWO 1.2 m and TRES 
observations, we proposed HAT-P-56 as a target for the 
K2 Campaign 0 through the Kepler Guest Observing 
Program. The observations are in Kepler Long Ca¬ 
dence mode (^ 30 min exposures) with a stamp size of 
27x27 pixels on Kepler CCD Module 10, Channel 29, and 
were carried out between BJD 2456728.5282 and BJD 
2456805.1883 (UT 2014 March 8 to UT 2014 May 27). 
There are two data gaps during the observation, from 
BJD 2456732.4309 to BJD 2456735.6386, and from BJD 
2456744.1180 to BJD 2456767.5941. Module 10 is one 
of the outermost modules on the Kepler spacecraft, and 
as a result, the target drifted for a significant fraction of 
a pixel during each 30 min exposure, leading to an elon¬ 
gated PSF. The bright neighbor 49" away from the target 
is partially observed in the target “postage stamps” (see 
Figure |5|). 

The data were reduced using HATNet’s reduc tion 
pipelin e, certain aspects of which are described in iPall 
(1200^ . After source extraction, we measured the flux 
of the star in a series of circular apertures, and the sky 
background in circular annuli. The sky background is 
determined by taking the median of all the pixels with 
iterative outlier rejection in the annuli to exclude the in¬ 
fluence from the bright neighbor star. The best aperture 
was of 3.3 pixels radius, based on the minimum r.m.s. we 
achieve (shown in Figure |S] together with the background 
annulus we used). 

To correct the photometry for variations due to the 
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TABLE 1 

Differential photometry of HAT-P-56. 


BJD^ 

(2,400,000-1-) 

Phot*’ 

o-phot Phot(orig)'’ 

Filter 

Instrument 

55954.98325 

0.00245 

0.00430 

r 

HATNet 

55971.72852 

0.00273 

0.00463 

r 

HATNet 

55860.09538 

0.00383 

0.00396 

r 

HATNet 

55957.77449 

0.00368 

0.00464 

r 

HATNet 

55929.86619 

0.00725 

0.00426 

r 

HATNet 

55901.95813 

-0.00271 

0.00432 

r 

HATNet 

55888.00414 

-0.01589 

0.00408 

r 

HATNet 

55904.75041 

0.00110 

0.00603 

r 

HATNet 

55929.86812 

0.00151 

0.00373 

r 

HATNet 

55901.95982 

0.00766 

0.00400 

r 

HATNet 


Note. — This table is available in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding its form and content. 
Data is also available at hatnet. org 

^ Barycentric Julian Date calculated directly from UTC, without correction for 
leap seconds for HATNet and KeplerCam. 

For HATNet and KeplerCam this is in units of magnitudes, for K2 it is in 
relative flux. The out-of-transit level has been subtracted. These values have 
been corrected for trends simultaneously with the transit fit for the follow-up 
data. For HATNet trends were filtered before fitting for the transit. 

^ Raw photometry values after correction using comparison stars, but without 
additional trend-filtering. For KeplerCam this is in magnitudes, for K2 it is in 
relative flux. We do not report this value for HATNet. 


TABLE 2 

Relative radial velocities, and bisector span measurements of 
HAT-P-56. 


BJD^’^ 

(2,456,900-1-) 

RV*" 

(m s~^) 

(ms-1) 

BS 

(m s“^) 

O-RS 
(m s~^) 

Phase 

Instrument 

32.99501 

68 

114 

48.70 

82.40 

0.937 

TRES 

34.00383 

-376 

152 

150.10 

57.10 

0.299 

TRES 

34.99731 

0 

54 

64.40 

58.30 

0.655 

TRES 

36.00811 

-310 

107 

2.40 

45.50 

0.017 

TRES 

42.99104 

-170 

70 

67.10 

50.00 

0.519 

TRES 

43.98771 

-268 

146 

-161.50 

56.80 

0.876 

TRES 

44.97264 

-558 

119 

-57.30 

70.80 

0.229 

TRES 

46.00112 

-113 

102 

-63.80 

74.60 

0.598 

TRES 

58.97630 

-493 

99 

-65.60 

44.80 

0.247 

TRES 

60.02365 

-221 

92 

-32.30 

65.70 

0.622 

TRES 

60.97679 

-150 

114 

49.20 

67.20 

0.964 

TRES 

61.96910 

-544 

78 

-111.70 

111.00 

0.319 

TRES 

65.96895 

-80 

87 

-43.50 

52.20 

0.753 

TRES 

69.93743 

-649 

95 

-2.40 

69.50 

0.174 

TRES 

70.89553 

-419 

153 

16.40 

58.90 

0.518 

TRES 

71.85396 

-187 

105 

85.90 

80.70 

0.861 

TRES 

72.94426 

-609 

54 

36.40 

34.10 

0.252 

TRES 

77.04183 

-331 

96 

-52.80 

47.40 

0.720 

TRES 


^ Barycentric Julian Date calculated directly from UTC, without correction for 
leap seconds. 

^ The zero-point of these velocities is arbitrary. An overall offset 7rei fitted to 
these velocities in Section |3] has not been subtracted. 

motion of the spacecraft, we performed an External Pa¬ 
rameter Decorrelation (EPD) on the extracted raw light 
curve. Before detrending, data points obtained during 
thrust fires were rejecte d following a similar methodol - 
ogy to that described in Vanderburg fc JohnsonI (|2014H . 

We also discarded the first segment of data (from BJD 
2456728.752975 to BJD 2456771.394346) obtained before 
the first safe mode, and the first transit right after the 
second safe mode while the space craft was still adjusting 
its pointing. We describe the systematics as a function 
of the X, y centroid of the star, the background flux bg, 


and the uncertainty of the background flux Cbg: 

f{m) = Co + Cl sin(27rai) -I- C2 cos(27rx) 

-I- C3 sin(27r y) + C4 cos(27r y) 

+ C 5 sin( 47 r a;) + cq cos( 47 r x) 

+ C 7 sin( 47 r y) + cs cos( 47 r y) 

-\- c^hg + cioCbg. 

To preserve the signal of the transit, we modify the 
original EPD algorithm to only fit with the out-of-transit 
part of the light curve. The light curve is then detrended 
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TABLE 3 

Stellar Atmospheric Parameters for HAT-P-56 Measured with SPC. 


BJD^ 

(2,456,900+) 

Teff.l 

(K) 

log 51 *' 
(c.g.s.) 

[m/H]i *> 

Dsinieff,! 
(km s“*) 

Teff.2 " 
(K) 

[m/H ]2 " 

'lisin2eff,2 
(km s“^) 

S/N " 

32.99501 

6406 

4.01 

-0.140 

40.06 

6545 

-0.052 

39.84 

55.7 

34.00383 

6500 

4.09 

-0.029 

40.72 

6592 

+0.034 

40.67 

56.2 

34.99731 

6500 

4.15 

-0.174 

40.19 

6538 

-0.150 

40.10 

102.7 

36.00811 

6588 

4.27 

-0.103 

39.83 

6576 

-0.116 

39.87 

84.0 

42.99104 

6528 

4.17 

-0.100 

39.62 

6572 

-0.079 

39.54 

92.9 

43.98771 

6527 

4.21 

-0.101 

40.20 

6553 

-0.099 

40.18 

83.8 

44.97264 

6533 

4.23 

-0.104 

39.86 

6531 

-0.100 

39.89 

80.1 

46.00112 

6551 

4.21 

-0.095 

39.75 

6569 

-0.091 

39.74 

101.6 

58.97630'* 





6555 

-0.071 

40.38 

91.9 

60.02365'* 





6598 

-0.057 

40.15 

74.9 

60.97679'* 





6534 

-0.109 

39.59 

90.6 

61.96910'* 





6544 

+0.062 

41.80 

49.3 

65.96895'* 





6576 

-0.082 

40.03 

113.2 

69.93743'* 





6603 

-0.077 

40.02 

91.9 

70.89553'* 





6612 

-0.083 

39.80 

78.6 

71.85396'* 





6540 

-0.093 

39.88 

97.6 

72.94426'* 





6592 

-0.123 

39.59 

96.5 

77.04183'* 





6561 

-0.099 

40.03 

78.4 


^ Barycentric Julian Date calculated directly from UTC, without correction for leap seconds. 

^ Measurements from initial SPC iteration in which log^* was allowed to vary. 

^ Measurements from second SPC iteration in which we fixed log^* to 4.235. 

^ Observations were obtained after running our initial SPC analysis. They were included, however, in our 
fixed-log^* SPC analysis. 



Fig. 4.— Top panel: RV measurements from FLWO 1.5m/TRES 
shown as a function of orbital phase, along with our best-fit circular 
model (solid line; see Table [Sj, and our best-fit eccentric model 
(dashed line). Zero phase corresponds to the time of mid-transit. 
The center-of-mass velocity has been subtracted. Second panel: 
Velocity O — C residuals from the best fit circular orbit. Third 
panel: Bisector spans (BS). Note the different vertical scales of the 
panels. 


with cosine filters as described in iHuang et al.l (I2013D . 
We allowed a minimum period of 0.5 days for the filters, 
to preserve the transit signal. 

The resulting K2 light curve is shown in Figure |6l to¬ 
gether with our best-fit transit model. The residuals have 


per-point r.m.s. of 0.17mmag. In order to search for ad¬ 
ditional transit signals, we removed the signal of HAT- 
P-56b from the K2 light curve according to the best fit 
transit parameters, and then reprocessed the light curve 
with the same detrending steps. We ran BLS on the de¬ 
trended residual light curve. There is no evidence for 
the presence of other transits. We also conducted a grid 
search in phase to look for the secondary eclipse with 
the transit period fixed. We could not find a significant 
detection of the secondary eclipse. We put an 1-cr up¬ 
per limit of 65 ppm (0.065 mmag) on the depth of the 
secondary eclipse. This is estimated by computing the 
weighted average of the variance of points in-transit and 
the variance of points out-of-transit. The 1-cr upper limit 
on the depth of the secondary eclipse we obtained from 
the K2 light curve is much smaller than the predicted sec¬ 
ondary eclipse depth (0.7 mmag ) if a ssuming the system 
is a blended binary (see Section 

3. ANALYSIS 

3.1. Stellar Parameters 

The adopted stellar atmospheric parameters of HAT- 
P-56 are based on our SPC analysis of the 18 high 
S/N TRES spectra (Section 12.21 and Table [S]). SPC 
was applied individually to each spectrum, and we take 
the error-weighted mean of the individual atmospheric 
parameter measurements. The adopted uncertainties 
(0.1 dex for logg*, 50K for Teff*, O.OSOdex for [Fe/H] 
and 0.50 km for ?;sini) reflect our estimate of the 
systematic errors in this method based o n observations 
of spe ctroscopic standard stars. Following iSozzetti et al.l 
(|2007f) we combine the resulting Tgff* and [Fe/H] mea¬ 
surements with the p* measurements from our joint anal¬ 
ysis of the RV and light curve data fSection l3.2|) . and 
compare them to the Yonsei-Yale stellar evolution mod- 
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X [asec] 

Fig. 5.— A image stamp from K2 Campaign 0 containing HAT- 
P-56. The neighbor at 49^^away from HAT-P-56 is partially ob¬ 
served at the bottom of the stamp. The inner most circle indicates 
the optimal aperture we selected that gave the best out of transit 
light curve r.m.s.. The annulus between the outer two circles is 
used to calculate the background. 

els (lYi et al.ll2001t the comparison is shown in Figure [7]), 
to determine the physical stellar parameters (mass, ra¬ 
dius, age, luminosity, etc.). As is often the case, we find 
that the resulting logg* value from this modeling dif¬ 
fers significantly from that determined through our ini¬ 
tial SPC analysis, and we therefore carried out a sec¬ 
ond iteration of SPC, fixing logg* to the isochrone-based 
value. Repeating the joint RV-blight curve analysis and 
the stellar evolution look-up, we find that the log 5 * value 
had converged, and therefore did not carry out any fur¬ 
ther iterations. The final parameters that we adopt for 
HAT-P-56 are listed in Table S] together with identify¬ 
ing information and catalog photometry. We find that 
the star HAT-P-56 has a mass of 1.296 ± 0.036 Mq, a 
radius of 1.428 ± 0.030 i?©, an age of 2.01 ± 0.35 Gyr, 
and is at a redde ning-corrected distan ce of 310.5 ± 7.1 pc 
(where we use the lCardelli et al1ll989l extinction law with 
Rv = 3.1). 

3.2. Global Modeling of RVs and Light Curves 

We carried out a joint analysis of the TRES RVs 
and th e HATNet, Kepler Cam and K2 light curves fol- 
lo winglBakos et al] (1201011 with modifications described 
bv iHartman et all (1201211 . The RVs are modeled using 
a Keplerian orbit , whil e the light curves are fit with 
a IMandel fc Agoll (1200211 transit model, with qu a dratic 
lim b darkening coeffici e nts a dopted from IClaretl (|2004[1 
and iClaret fc BloemenI (1201111 . The HATNet light curve 
that we analyzed has been filtered via the EPD and TEA 
procedures before fitting the model, so we include a di¬ 
lution factor to account for possible over-filtering of the 
transits. For the K2 light curve, we integrate the model 
over the 30 min exposure time (this is done by evalu¬ 
ating the model flux ratio at four evenly spaced times 
within a 30 min bin, and taking the a verage). We used 
equation 8 of lKipping fc BakosI ((2011all for our K2 model 
to account for a possible occultation signature and/or a 
variation due to reflected light (in practice we find that 
the light curve is consistent with no such variations). 
The K2 light curve has been fitted simultaneously with 
the TFA algorithm using 47 template stars observed in 
the same channel. For the KeplerCam light curve, we 
use a simple model for instrumental trends (consisting of 


a quadratic function of time and linear functions in the 
point spread function (PSF) shape parameters), which 
we fit simultaneously with the physical model. We use a 
differe ntial evolution M arkov-chain Monte Carlo (DEM- 
CMC: Iter BraakI [200611 procedure to explore the fitness 
landscape and to determine the posterior parameter dis¬ 
tributions. 

The fit was performed both fixing the eccentricity to 
zero, and allowing it to vary. We find that the TRES 
RVs do prefer a slight eccentricity of e = 0.130 ± 0.058, 
but based on the Bayesian evidence (est imated from the 
Mark ov-chain following the method of IWeinberg et al.l 
EoH), we conclude that the difference in is not signifi¬ 
cant enough to justify the additional free parameters, and 
therefore we adopt the circular orbit model. The 95% 
confidence upper limit on the eccentricity is e < 0.246. 

Table [5] lists the adopted parameters for the planet 
HAT-P-56b. We find that this planet has a mass of 
2.18±0.25 Mj, a radius of 1.466±0.040i?j, and is orbiting 
its host star with a period of 2.7908327 ± 0.0000047days 
and an orbital separation of 0.04230±0.00039 AU. At this 
separation the planet would have an equilibrium temper¬ 
ature of 1840 ± 21K assuming zero albedo and complete 
redistribution of heat. 

3.3. Blend Analysis 

In order to rule out the possibility that HAT-P-56 is 
a blended stellar eclipsing binary system we carr ied out 
a blend analysis following IHartman et al.l (|2012ll . with 
a few modifications to properly handle the K2 light 
curve. These include: (1) integrating each simulated 
light curve model over the 30 min exposures and using 
the integrated m odel in calculating the difference from 
the observations (iKioping fc Bakosll2011all : (2) using the 
iClaret fc BloemenI (120111) limb darkening coefficients for 
the Kepler band-pass; (3) using a polynomial transfor¬ 
mation from griz to the Kepler Kp magnitude system to 
predict the relative fluxes of blended stars in the Kepler 
band-pass. 

We find that we can rule out a blended eclipsing binary 
scenario with greater than 5a confidence based solely on 
the photometry (including both HATNet and K2 pho¬ 
tometry). The K2 data is key in making this assessment. 
Although blend models exist which fit the primary tran¬ 
sit, such models predict a secondary eclipse with a depth 
of AKp = 0.7 mmag, which is ruled out by the K2 light 
curve (Figure [5]). Further evidence against a blend sce¬ 
nario is the lack of BS variations (the TRES BS mea¬ 
surements have an r.m.s. scatter of 77ms“^) and the 
significant RV variation in phase with the photometric 
ephemeris and consistent with a transiting planet. 

We also considered the possibility that HAT-P-56 is a 
transiting planet system with a fainter, unresolved, stel¬ 
lar companion. From the 2MASS catalog there are no 
known stars within 20" of HAT-P-56 with AK < 5 mag. 
The highest spatial resolution observations available are 
our KeplerCam observations which have a PSF FWHM 
of 3.5". Based on these images we can rule out a com¬ 
panion with Ai < 3.5 mag to within 5.4". Based on our 
blend modeling we find that models including a stellar 
companion with M > 0.67 Mq yield a higher value 
than a single star with a transiting planet, but, except 
for companions very close in mass to HAT-P-56b, which 
can be ruled out with > 3cr confidence, the difference in 
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Fig. 6 .— Phase-folded 30 min cadence K2 observations of HAT-P-56. The phase-folded light curve is shown in the bottom panel. In 
the top left we show a zoom-in on the transit, while in the top right we show a zoom-in on the expected phase of occultation assuming a 
circular orbit. The solid line shows our best-fit model to the observations, binned over the 30 min integrations. The dashed line, visible 
only in the top-left panel, shows the model without binning. 
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Fig. 7. — Comparison between the measured values ofT^ff* 
and p ir (filled circle), and the model isochrones from lYi et al.l 
1120011) . The best-fit values, and approximate Icr and 2(7 confidence 
ellipsoids are shown. The open triangle shows the values from our 
initial SPC iteration. The isochrones are shown for ages of 0.2 
to 3.0 Gyr, in 0.2 Gyr increments. 

between models with and without companions is not 
statistically significant. Higher spatial resolution imag¬ 


ing, and/or long term RV monitoring is needed to check 
for a binary star companion. If such a companion exists, 
the mass and radius of HAT-P-56b would both be larger 
than what we infer here. 

3.4. Out of transit variation of HAT-P-56 

Stellar variability with amplitude ^1 mmag are 
present in the K2 light curve (see Fi gure [9|. The Lomb- 
Scarg le normalized power spectrum (|Lomblll97fiI iScar in 
119821! of the star suggest that it is likely a pulsating 7 
Dor, with the primary period of ~ 1.644 ± 0.03 days. 
The second strongest period peak is at ~ 1.744 ± 0.023 
days. The effective temperature and the surface grav¬ 
ity of the star place it outside the low temperature 
boundary of the classical instability strip, but within the 
range of other 7 Dors disc overed by the Kepler Mission 
(|Uvtterhoeven et al.ll 2011 l !. 

We tried to constrain the orbital phase variation due 
to HAT-P-56b with the K2 light curve. To take into 
account the influence of the stellar variability, we simul¬ 
taneously fit for the amplitude of the six most dominant 
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Fig. 8 . — Comparison between the phased K2 observations of HAT-P-56 and the predicted light curve for the best-fit blended stellar 
eclipsing binary model (top curve and data in each panel) together with the residuals from the model (bottom of each panel). In the 
left panel we show the primary eclipse, while in the right we show the secondary eclipse. The simulated primary eclipse fits the observed 
transit, but a significant secondary eclipse is predicted, but not observed. Note the different vertical scales in the two panels. 
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Fig. 9. — The top panel shows a segment of the K2 light curve 
of HAT-P-56, demonstrating the out of transit variability pattern. 
The bottom panel shows the Lomb-Scargle normalized power spec¬ 
trum of the K2 light curve excluding the transits. 

Fourier modes detected by the Lomb-Scargle algorithm, 
together with the expected reflection, beaming, and ellip¬ 
soid variation effects with periods constrained due to the 
presence of the planet. The best fitted values and error 
bars are presented in Table El We report a 3 -ct detection 
of the reflection effect, with an amplitude of 21 ± 7 ppm, 
which is comparable to the theoreticall y estim ated value 
(27.4ppm, following iMazeh fc Faigleil (1201011 1. We do 
not detect the beaming effect. The ellipsoidal variation is 
detected with a 2-cr significance. However, we caution the 
reader, that we also obtain a high amplitude coefficient 
for the sin (27r/(Porb/2)t_,) term, which indicates the de¬ 
tected ellipsoidal variation is not in phase with the planet 
transit. This could be due to our poor understanding of 
the stellar variability, suggesting the modeled ellipsoidal 
variation amplitude may not be physically meaningful. 

4. DISCUSSION 

In this paper we presented the discovery and charac¬ 
terization of HAT-P-56b, an inflated massive hot Jupiter 
around a bright F star. 

The radius anomaly of hot Jupiters is on e of the oldest 
unsolved problems in the exop lanet field (|Baraffe et al.l 
l20irt [Siegel fc Burrowsll2012l ). Observationally, only a 
few very inflated planets {R > 1.35i?j) are observed to be 
above the mass of 2Mj. It has been empirically demon¬ 
strated by variou s authors that more massive planets are 
harder to inflate. lEnoch et al.l (120121) derived the depen¬ 


dence of giant planet radius in three mass regimes, and 
found that the radii of high mass planets (M > 2Mj) are 
less sensitive to the equilibrium temperature of the plan¬ 
ets, a s compared to Jupiter mass planets. I Weiss et aH 
(|2013ll fitted for the fundamental plane of giant planets 
using existing data. They found that the radius of irra¬ 
diated hot Jupiter s inver sely correlates with the planet 
mass. I Zhou et all (j2014l l investigated the mass depen¬ 
dence of planet radius on the equilibrium temperature, 
and also found that as the planet mass increase above 
IMj, the influence on radius from irradiation decreases. 

In Figure [TO] we show the location of HAT-P-56b in 
the mass radius diagram of hot Jupiters. Taken at face 
value, it is the most inflated hot Jupiter with a mass be¬ 
tween 1.5Mj—4Mj. We show the position of HAT-P-56b 
relativ e to the fundamental plane fitted bv iWeiss et al.l 
(1201311 in Figureim where it can be seen that HAT-P-56b 
is the farthest away from the fundamental plane when 
compared to the other massive planets (green points in 
the figure). The fundamental plane of irradiated hot 
Jupiter s (M>150Mm, or 0 .47M.t), is expressed as the fol¬ 
lowing (|Weiss et al.ll2013ll : 


= 2.45(|)£)-»“( 

R® -'M'e 


F 




0.094 


ergs ^cm ^ 


The unusually large radius of HAT-P-56b, given its mass, 
makes it an important data point in the hot Jupiter pop¬ 
ulation. 

However, it is noteworthy that the radius of HAT-P- 
56b may have a large uncertainty due to the large im¬ 
pact parameter (6 = 0.8725^0^0060) indicating a nearly 
grazing transit. The grazing transit geometry makes 
the transit depth highly dependent on the limb dark¬ 
ening parameters of the star. The radius uncertainties 
are largely reduced thanks to the high precision K2 light 
curve. Another way to better constrain the radius of the 
planet would be to observe transits in the near infrared 
where limb darkening is negligible. The only other plan¬ 
ets known to be on close-to grazing orbits are TrES-2b, 
HAT-P-27b, WASP-34b . WASP-67b, and Kepler-447b 
(lO’Donovan et all I2006I iBekv et ahl I20I1 : FHellier et al.l 


l20I2t iMancini et al.ll2014l iLillo-Box et ah' l20I5f) . While 

the grazing configuration increases the uncertainty o f 
the planetary radius measurement, as iRibas et al.1 (|2008[1 
pointed out, such a near grazing transit has the advan- 
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TABLE 4 

Stellar Parameters for HAT-P-56 


Parameter 

Value 

Source 

Identifying Information 

R.A. (h:m:s) 

06*'43™23.52s 

2MASS 

Dec. (d:m:s) 

-|-27°15'08.2" 

2MASS 

GSC ID 

GSC 1901-00976 

GSC 

2MASS ID 

2MASS 06432353-1-2715082 

2MASS 

EPIC ID 

202126852 

EPIC 

Spectroscopic properties 

Teff* (K). 

6566 ± 50 

SPC ^ 

[m/H]. 

-0.077 ±0.080 

SPC 

V sinf (kms ^) ... 

40.06 ±0.50 

SPC 

7 rv (kms-1). 

35.11±0.1 

TRES 

Photometric properties 

B (mag). 

11.287 ±0.040 

APASS 

V (mag). 

10.908 ± 0.036 

TASS Mark IV 

(mag). 

9.919 ±0.093 

TASS Mark IV 

Kep (mag). 

10.9 

EPIC 

g (mag). 

11.040 ±0.030 

APASS 

r (mag). 

10.762 ±0.080 

APASS 

* (mag). 

10.709 ±0.060 

APASS 

J (mag). 

10.068 ± 0.021 

2MASS 

H (mag). 

9.881 ±0.017 

2MASS 

Ks (mag). 

9.830 ±0.016 

2MASS 

Derived properties 

M* (Mo). 

1.296 ±0.036 

Isochrones±p*±SPC 

. 

1.428 ±0.030 

Isochrones+p*+SPC 

log a* (cgs). 

4.240 ±0.015 

Isochrones+p*+SPC 

L. (Do). 

3.39 ±0.19 

Isochrones+p*+SPC 

Mv (mag). 

3.411 ±0.067 

isochrones+p*+SP(J 

Mk (mag.ESO) 

2.398 ±0.048 

Isochrones+p*+SPC 

Age (Gyr). 

2.01 ±0.35 

Isochrones+p*+SPC 

Av (mag) ==. 

0.00801“;°5» 

Isochrones+p*+SPC 

Distance (pc). 

310.5 ± 7.1 

isochrones+p*+SP(J 

P*(g cm-3). 

0.627 ±0.033 

LC 


^ SPC = “Stellar Parameter Classification” method based on cross-correlating high- 
resolution spectra against synthetic templates iBuchhave et al.ll2012D . These parame¬ 
ters rely primarily on SPC, but have a small dependence also on the iterative analysis 
incorporating the isochrone search and global modeling of the data, as described in 
the text. 

^ Isochrones-|-p*-|-SPC = Based on the isochrones er et al.ll200il> . the stellar 
density used as a luminosity indicator, and the SPC results. 

^ Total V band extinction to the star determined by comparing the catalog 
broad-band photometry listed in the table to the expected magnitudes from the 
Isochrones-|-p*-|-SPC model for the star. We use the lCardelli et al.l (119891) extinction 
law. 

^ The stellar density is determined primarily from fitting the transit light curve. 
There is also a slight dependency on the SPC atmospheric parameters and stellar 
evolution models through the adopted limb darkening coefficients. 


tage of its depth and duration being more sensitive to the 
presence of other planetary companions on inclined or¬ 
bits. This makes HAT-P-56b a promising target among 
the hot Jupiters for detecting transit timing and duration 
variations. 

Given the brightness of the host star (Lnag = 10.9), 
HAT-P-56 is a system of particular interest for measur¬ 
ing the spin-orbit obliquity angle. The estimated RM 
effe ct amplitude A Vr is ^ 140ms“^ following Equation 
5 of IGaudi fc WinnI (120071 1 . HAT-P-56 is a relatively fast 
rotating F star w ith effective teinperat ure (6566 ± 50 K) 
slightly above the (Albrecht et ahl (j2012ll division for tidal 
alignment of hot Jupiters. Measurement of the spin- 
orbit obliquity will contribute an important data point 
to the obliquity distribution statistics. The rotation of 
HAT-P-56 is similar to GoRoT-11, and rapid enough that 


the distortion in the spectral line profiles during transit 
may be reso lved, enabling transit Doppler tomography of 
the s ystem (IGandolfi et al.l[20T^ (Collier Cameron et al.1 
(201011 . Similar to CoRoT-11, the rotation period of HAT- 
P-56 is fast er than the orbit period of HAT-P-56b, which 
is unusual ((Walkowicz fc Basril(2013ll . This may either 
suggest that the tidal interaction between the star and 
planet is weak, or the planet might be pushed out by the 
tides of the star from a closer orbit. 
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by NASA grants NNG04GN74G and NNX13AJ15G. 
Follow-up of HATNet targets has been partially sup¬ 
ported through NSF grant AST-1108686. G.A.B., Z.C. 
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TABLE 5 

Parameters for the transiting planet HAT-P-56b. 


Parameter 

Value ^ 

Light curve parameters 


P (days) . 

2.7908327 ±0.0000047 

Tc (BJD) b . 

2456553.61645 ± 0.00042 

Ti 4 (days) . 

0.09463 ± 0.00093 

Ti 2 = T34 (days) . 

0.0336 ± 0.0021 

a/Ri, . 

6.37 ±0.11 

C/R* . 

29.37 ±0.21 

Rp/R* . 

0.10540 ±0.00086 

. 

U. /D10_Q QJ^Q4 

b = a cos i / R* . 

n q79c:-I-0.0044 

U.O/ZO_Q QQQQ 

* (deg) . 

Limb-darkening coefficients ^ 

82.13 ± 0.18 

Cl,2 (linear term) . 

0.1716 

C2,2 (quadratic term) ... 

0.3692 

ci,r . 

0.2396 

C2, r . 

0.3825 

Cl, Kep . 

0.3239 

C2, Kep . 

0.3235 

RV parameters 


K (ms-l) . 

262 ± 30 

e ^ . 

< 0.246 

RV jitter (ms“^) ^ . 

< 15 

Planetary parameters 


Mp (Mj) . 

2.18 ±0.25 

Rp (Rj) . 

1.466 ± 0.040 

C{Mp,Rp) g . 

0.08 

Pv (gcm-3) . 

0.86 ±0.12 

loggp (cgs) . 

3.402 ± 0.055 

a (AU) . 

0.04230 ± 0.00039 

Req (K) h . 

1840 ± 21 

0 ‘ . 

0.096 ± 0.011 

{Fj) (lO^ergs-l cm~^) ' 

2.59 ±0.12 


^ The adopted parameters assume a circular orbit. Based 
on the Bayesian evidence ratio we find that this model is 
strongly preferred over a model in which the eccentricity is 
allowed to vary in the fit. For each parameter we give the 
median value and 68.3% (l*^) confidence intervals from the 
posterior distribution. 

° Reported times are in Barycentric Julian Date calculated 
directly from UTC, without correction for leap seconds. Tc'. 
Reference epoch of mid transit that minimizes the correlation 
with the orbital period. T 14 : total transit duration, time 
between first to last contact; T 12 = T 34 : ingress/egress time, 
time between first and second, or third and fourth contact. 

^ Reciprocal of the half duration of the transit used as a 
jump parameter in our MCMC analysis in place of a/R*. It 
is related to a/R* by the expression Q/R* — a/R*(27r(l + 
e sin a;))/(P'\/l — b'^y/1 — e^) JBakos et al.ir2010n . 

^ Values for a quadratic law, adopted from the tabulations 
bv IClar^ j2004D : IClaret et al.| l |2013l) according to the spec¬ 
troscopic (SPC) parameters listed in Table 
® The 95% confidence upper-limit on the eccentricity from a 
model in which the eccentricity is allowed to vary in the fit. 
^ Error term, either astrophysical or instrumental in origin, 
added in quadrature to the formal RV errors. This term is 
varied in the fit assuming a prior inversely proportional to 
the jitter. In this case we find a preferred value of 0 for the 
jitter, and list the 95% confidence upper limit. 

® Correlation coefficient between the planetary mass Mp 
and radius Rp determined from the parameter posterior dis¬ 
tribution via C(Mp,Rp) —< (Mp— < Mp >)(Rp— < Rp > 
) > /{<yMp<yRp) > where < • > is the expectation value 
operator, and (7x is the standard deviation of parameter x. 

Planet equilibrium temperature averaged over the orbit, 
calculated assuming a Bond albedo of zero, and that flux is 
reradiated from the full planet surface. 

’ The Safronov number is given by © — -^(V esc/Vorb)^ = 
(a/Rp)(Mp/M* ) fsee lHansen Sz Barman||2007D . 

J Incoming flux per unit surface area, averaged over the 
orbit. 
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TABLE 6 __ 

Derived and E xpecte d Orbit Phase Variation UMazeh fc Faigl^ 
1201011 FOR HAT-P-56 light curve 


Coefficient 

Derived Value 
(ppm) 

Expected Value ® 
(ppm)=‘ 

Effect 

aic 

-21I7 

-27.4 ±0.6 

Reflection 

ais 

istl 

3.5 ±0.4 

Beaming 

a2C 

-17+7 

-6.5 ±0.8 

Ellipsoidal 

a2S 

45^7 

- 

- 


^ The expected values are computed following lMazeh fc Faigleil J201C1I1 . 
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Fig. 10.— Mass-radius diagram for HAT-P-56b (big square) 
compared to the full sample of confirmed hot Jupiters (circles) 
from http://exoplanets.org The points are color coded with their 
equilibrium temperature (cooler planets have darker colors). The 
err or bars indicate one s igma uncertainties. The lines represents 
the IFortnev et al.l 1120071') planet radius model without heavy ele¬ 
ments cores. The black lines are for planets at 0.045 AU from a sun 
like star (solid: IGyr, dashed:3 Gyr). The red lines are for planets 
at 0.02 AU from a sun like star (solid:lGyr, dashed:3 Gyr). 
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Predicted Radius (Rj) 

Fig. 11.— Observed radius of hot J upiters versus the radius 
predicted by the empirical relation of (I Weiss et al.ll20T^ . The 
green filled circles represent planets with mass 1.5Mj < Mp < 
3.0Mj; the light grey open triangles represent planets with mass 
0.47 < Mp < 1.5Mj. The red square is HAT-P-56b. 
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